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Abstract. We provide a complete description of the critical threshold phenomena for 
the two-dimensional localized Euler-Poisson equations, introduced by the authors in ]l7p . 
Here, the questions of global regularity vs. finite-time breakdown for the 2D Restricted 
Euler-Poisson solutions are classified in terms of precise explicit formulae, describing a 
remarkable variety of critical threshold surfaces of initial configurations. In particular, it 
is shown that the 2D critical thresholds depend on the relative size of three quantities: 
the initial density, the initial divergence as well as the initial spectral gap, that is, the 
difference between the two eigenvalues of the 2x2 initial velocity gradient. 
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1. Introduction and statement of main results 

We are concerned with the critical threshold phenomena in multi-dimensional Euler- 
Poisson equations. In this paper we consider a localized version of the following two- 
dimensional Euler-Poisson equations 

)2 



(1.1) 

(1.2) 
(1.3) 



d t p + V • (pU) = 0, xgIR", tGlFV 
dt(pU) + V-(pU®U) = -kpV0, 
- A(f) = p - c, x e IR 2 , 



which are the usual statements of the conservation of mass, Newton's second law, and 
the Poisson equation defining, say, the electric field in terms of the charge. Here k > 
is a scaled physical constant, which signifies the property of the underlying repulsive 
forcing (avoiding the case of an attractive force with k < 0), and c denotes the constant 
"background" state. The unknowns are the local density p = p(x,t), the velocity field 
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U = (u,v)(x,t), and the potential <fi = <f>(x,t). It follows that as long as the solution 
remains smooth, the velocity U solves a forced transport equation 

(1.4) d t U + U ■ VU = F, F = -kV(j) 

with being governed by Poisson's equation (|1.3|) . 

This hyperbolic-elliptic coupled system (|l.lp-(|IT3"D describes the dynamic behavior of 
many important physical flows including the charge transport [22|, plasma with collision 
i~3| , cosmological waves || and the expansion of the cold ions [12]. Let us mention that the 



Euler-Poisson equations could also be realized as the semi-classical limit of Schrodinger- 
Poisson equation and are found in the 'cross-section' of Vlasov-Poisson equations. These 
relations have been the subject of a considerable amount of work in recent years, and we 
refer to ]1(J, || and references therein for further details. 

To put our study in a proper perspective we recall that there has been a considerable 
amount of literature available on the global behavior of Euler-Poisson and related prob- 
lems, from local existence in the small if s -neighborhood of a steady state flTj|, 23, [J to 
global existence of weak solutions with geometrical symmetry ||, for the two-carrier types 



in one dimension [29|, the relaxation limit for the weak entropy solution, consult [21] for 



isentropic case, and |]T4j for isothermal case. 

For the question of global behavior of strong solutions, however, the choice of the initial 
data and/or damping forces is decisive. The non-existence results in the case of attractive 
forces, k < 0, have been obtained by Makino-Perthame [20|, and for repulsive forces by 
Perthame [23 1. For the study on the singularity formation in the model with diffusion 
and relaxation, consult |HJ. In all these cases, the finite life span is due to a global 
condition of large enough initial (generalized) energy, staying outside a critical threshold 
ball. Using the characteristic-based method, Engelberg [0 gave local conditions for the 
finite-time loss of smoothness of solutions in Euler-Poisson equations. Global existence 
due to damping relaxation and with non-zero background can be found in 
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For the model without damping relaxation the global existence was obtained by Guo [|Tl 



assuming the flow is irrotational. His result applies to if 2 -small neighborhood of constant 
state. 

When dealing with the questions of time regularity for Euler-Poisson equations without 
damping, one encounters several limitations with the classical stability analysis. Among 
others issues, we mention that 

(i) the stability analysis does not tell us how large perturbations are allowed before 
losing stability - indeed, the smallness of the initial perturbation is essential to make the 
energy method work, e.g. 
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(ii) the steady solution may be only conditionally stable due to the weak dissipation in 
the system, say in the ID Euler-Poisson equations ||. 

In order to address these difficulties we advocated, in H, a new notion of critical 
threshold (CT), which describes the conditional stability of the ID Euler-Poisson equa- 
tions, where the answer to the question of global vs local existence depends on whether 
the initial configuration crosses an intrinsic, 0(1) critical threshold. Little or no atten- 
tion has been paid to this remarkable phenomena, and our goal is to bridge the gap of 
previous studies on the behavior in Euler-Poisson solutions, a gap between the regularity 
of Euler-Poisson solutions in the small and their finite-time breakdown in the large. The 
critical threshold (CT) in the ID Euler-Poisson system was completely characterized in 



Critical Thresholds in restricted Euler-Poisson equations 



3 



terms of the relative size of the initial velocity slope and the initial density. Moving to 
the multi-D setup, one has first to identify what are the proper quantities which govern 
the critical threshold phenomena . In |l7j we have shown that these quantities depend in 
an essential manner on the eigenvalues of the gradient velocity matrix, Vtt. In order to 
trace the evolution of M := VC/, we differentiate Q1.4Q , obtaining formally 

(1.5) d t M + U ■ VM + M 2 = -jfe(V ® V)0 = kR[p - c], 

where R[ } is the 2x2 Risez matrix operator, defined as 



R[f] =: V ® VA- 1 ^] = T' 1 [^f(0 



j,k=l,2 



The above system is complemented by its coupling with the density p which is governed 

(1.6) d t p + U ■ Vp + ptrM = 0. 

Passing to the Lagrangian coordinates, that is, using the change of variables a i— > x(a, t) 
with x(a, t) solving 

— = U(x,t), x(a, 0) = a, 
Euler-Poisson equations are recast into the coupled system 

(1.7) ^M + M 2 = kR[p-c], 

(1.8) ^p + ptrM = 0, 

at 

with d/dt standing for the usual material derivative, d t + U ■ V. It is the global forcing, 
kR[p — c], which presents the main obstacle to study the CT phenomena of the multi-D 
Euler-Poisson setting. 

In this work we focus on the Restricted Euler-Poisson (REP) system introduced in 
||17|| , which is obtained from ( |1 . 7| ) by restricting attention to the local isotropic trace, 
|(p — c)/2x2, of the global coupling term kR[p — c], namely 

(1.9) ^M + M 2 = |(p-c)-/ 2x2 , 

(1.10) jp + ptrM = 0. 

We are concerned with the initial value REP problem ( |1.9|) ,( |ITT0D , subject to initial data 

(M,p)(-,0) = (M ,po). 

We note in passing that the REP system is to the full Euler-Poisson equations what 
the Restricted Euler model is to the full Euler equations, consult |26|, |3], [[], £|, |T7|]. The 
existence of a critical threshold phenomena associated with this 2D REP model with zero 
background, c = 0, was first identified by us |]T7J. The current paper provides a precise 
description of the critical threshold for the 2D REP system ( |1.9| ), QLlUp , with both zero 
or non-zero background charge. In particular, we use the so called Spectral Dynamics 
Lemma, [17, Lemma 3.1] to obtain remarkable explicit formulae for the critical threshold 
surfaces summarized in the main Theorems ITTTl and |L2| below. 
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To state our main results, we introduce two quantities with which we characterize the 
behavior of the velocity gradient tensor M. These are the trace, d := trM, and we note 
that in case M coincides with VU then d stand for the divergence, d = u x + v y ; and the 
nonlinear quantity T := (trM) 2 — AdetM, which serves as an index for the spectral gap. 
Indeed, if Aj, % — 1, 2 are the eigenvalues of M, then 

\i = \[d-Vf], \ 2 = ^[d + Vf] 

and hence T is nothing but the square of the spectral gap T = (X 2 — Xi) 2 . We note that 
when M coincides with VZ7, then T = (u x — v y ) 2 + 4u y v x , and the role of this spectral 
gap was first identified in the context of the 2D Eikonal equation in [17, lemma 5.2]. 



We observe that if T < then the spectral gap is purely imaginary, and otherwise the 
spectral gap is real. 

Theorem 1.1. [2D REP with zero background] . Consider the 2D repulsive REP system 
( \1.9j )- ([T7n] ) with k > and with zero background c = 0. Then the solution of 2D REP 
remains smooth for all time if and only if the initial data (pq,Mq) lies in one of the 
following two regions, (p , d , T ) G Si U S%: 
(i) Either 



T < and 



d>0 if p = 

d arbitrary if p > 



(po,do,r ) g Si, Si := |0,d,r) 

(ii) or 

Oo, d , T ) G S 2 , S 2 := d,T) p > 0, T > 0, and d>g[p,T)} 



where 



g{p,T) :— sgn(r — 2kp) \ T — 2kp + 2kpln 



2kp 



Theorem 1.2. [2D REP with nonzero background]. Consider the 2D repulsive REP sys- 
tem ( \l.Q )- flTM\ ) with k > and with non-zero background c > 0. Then the solution of 
2D REP remains smooth for all time if and only if the initial data (po, M ) lies in one of 
the following three regions, (p , d , T ) G Si U S 2 U S 3 : 
(i) Either 



T < and 



d>0 



ifp = 



(po, d , T ) G S±, Si := <{ (p, d, T) 

(ii) or 

(po, do, T )eS 2 , S 2 := | (p, d,T) < T < and | ^>gl(p,T) ' if T > 2k(p 

where 



k 



d arbitrary if p > 



\d\< gi (p,T) if T<2k(p-c) 

c) 



9i(p,T) 



T-2k 



+ Vp 2 - 2ck~ x T + pi 



71 



p- y/p 2 - 2ck- 1 T 



2c 
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(Hi) or 



f />,,. f/.i. i n ) e -s' :i . ,s' :! : . \ f (/. n r = ^- P 2 , d = g 2 (p, r), p>o} . 



where 



k n , , /2c 



0a(/O = ^i(p,r)| r= ^ p2 := J -2ck + —p 2 + 2kpln [ — 



P 



Several remarks are in order. 

1. The above results show that the global smooth solution is ensured if the initial 
velocity gradient has complex eigenvalues, which applies, for example, for a class of initial 
configurations with sufficiently large vorticity \uq v — Vq x \ » 1. With other initial config- 
urations, however, the finite time breakdown of solutions may - and actually does occur, 
unless the initial divergence is above a critical threshold, expressed in terms of the initial 
density and initial spectral gap. Hence, global regularity depends on whether the initial 
configuration crosses an intrinsic, 0(1) critical threshold. 

2. The critical threshold in the ID Euler-Poisson equations depends on the relative 
size of the initial velocity slope and the initial density, consult ||. In contrast to the ID 
scenario, the critical threshold presented here depends on three initial quantities: density 
p , divergence V • U and initial spectral gap T = (u 0x — v 0y ) 2 + 4uo y v 0x . 



3. Theorem |1.1| tells us that the size of initial sub-critical range which gives rise to 
regular solution is decreasing as the initial ratio To/po is increasing. In particular when 
this ratio is larger than 2k, then the initial divergence must stay above a positive critical 
threshold to avoid the finite time breakdown. 



4. From Theorem |L2] we see that the initial critical range which guarantees global reg- 
ularity shrinks as the initial ratio r /po is increasing in (— oo, ^-). Finite time breakdown 
must occur when this ratio is larger than j-. 

5. The limit c j is a sort of a singular limit and hence one cannot recover Theorem 



simply by passing to the limit c — ► in Theorem |1.2| . □ 
It is well known that a finite time breakdown is a generic phenomena for nonlinear hy- 
perbolic convection equations, which are realized by the formation of shock discontinuities. 
In the context of Euler-Poisson equations, however, there is a delicate balance between the 
forcing mechanism (governed by Poisson equation), and the nonlinear focusing (governed 
by Newton's second law), which supports a critical threshold phenomena. 

In this paper we show how the persistence of the global features of the solutions for 
REP hinges on a delicate balance between the nonlinear convection and the localized 
forcing mechanism dictated by the Poisson equation. Here we use these restricted models 
to demonstrate the ubiquity of critical thresholds in the solutions of some of the equa- 
tions of mathematical physics. This remarkable CT phenomena has been found in other 
contexts, say convolution model for nonlinear conservation laws [16j , a nonlocal model in 
the nonlinear wave propagation P5|, etc. 



We note in passing that in this paper we restrict ourselves to Euler-Poisson equations 
with localized forcing. The presence of the global forcing allows for additional balance, 
and we hope to explore the critical threshold phenomena for the general model with global 
forcing in a future work. 
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We now conclude this section by outlining the rest of the paper. In Section 2 we study 
the critical threshold for the REP with zero background. The key observation is that 
the spectral gap is conserved along particle path. With this property we will be able to 
reduce the full dynamics on the 2D manifold parameterized by this initial spectral gap. In 
Section 3 we discuss the critical threshold for the REP with non-zero background, where 
the CT arguments become considerably more involved. We treat the different cases which 
are indexed by the initial spectral gap. 

2. 2D Restricted EP with zero background 

In this section we prove the existence of the critical threshold of the 2D restricted EP 
with zero background (c = 0) 



d , , . , k 

"'2x2, 



(2.1) jM + M 2 = |p/ 2xz . 

(2.2) 4p + pt rM = 0. 

dt 

This system with initial data (p , M ) is well-posed in the usual H s Sobolev spaces for 
a short time. The global regularity follows from the standard boot-strap argument once 
a priori estimate on ||M(-)||£°° is obtained. First we show that for the 2D restricted EP 
(|2.1D-(PT2"D, the velocity gradient tensor is completely controlled by the divergence d and 
the density p. 

Lemma 2.1. Let M be the solution of the 2D restricted EP, then the bounded of M 
depends on the boundedness oftrM and p, namely, there exist constant, Const = Constx 
such that 

||M(-,i)IU-[o,T] < Const T . || (trM,p) || L ~ [0iT] . 

Proof. For the 2-D case the velocity gradient tensor is completely governed by p := My — 
M 22 , q := M12 + M 21 , u = M l2 - M 2 i and d = M u + M 22 . From the M-equation 



d ( M n Afia \ ( Ml + M 2l M 12 dM l2 \ k 



'2x2, 



dt V M 2 i M 22 J + \ dM 2l M 21 M 12 + Ml 2 ) 2 pI ' 

one can obtain 

d J n 

— p + p d = 0, 

Jt q + qd = ^ 
d 

—us + cud = 0, 
dt 

which when combined with the mass equation 

d 

dt P + pd = ° 

gives 

(p,q,u) = (p ,qo,Uo)Po 1 P- 
This shows that | My 1 2,00 are bounded in terms of |g?|z,oo and \p\l^ as asserted. □ 
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This lemma tells us that to show the global regularity it suffices to control the divergence 
d and the density p. Let Aj,z = 1,2, be the eigenvalues of the velocity gradient tensor, 
then d = Ai + A 2 and the continuity equation Q1.10Q reads 



(2.3) ^ p + p ( Al + A 2 )=0. 



The Spectral Dynamics Lemma fT7 , Lemma 3.1] tells us that the velocity gradient equa- 
tion fln]D yields 

(2.4) ^ 1 + A-|p, 

(2.5) & + % = 



Following [0 we consider the difference of the last two equations which gives for 77 : = 
A2 — Ax, 

- v + rj(X 1 + X 2 ) = 0. 
This combined with the mass equation (|2.3|) and trM = Ai + A 2 yields 

dt \pj p po[a) 

Set (3 := ?7g(a)/pg(a) as a moving parameter with the initial position a G IR 2 , one then 
obtains a closed system for p and d 

(2.6) p' + pd = 0, ':=^, 

(2.7) d + = kp. 



The first is the mass equation; the second is a re-statement of summing (273),(E73), d' + 
(d 2 + r] 2 )/2 = kp with r] 2 = (3p 2 . 

We shall study the dynamics of (p, d) parameterized by (3. It is easy to see that if 
the initial eigenvalues are complex, then the eigenvalues remain complex as time evolves. 
From 

= % r = (A 2 (0) - A x (0)) 2 
Po 

we see that we need to distinguish between two cases, namely, (3 < where the initial 
spectral gap is complex, and (3 > where the initial spectral gap is real. 



2.1. Complex Spectral Gap. We first study the case (3 < when the initial eigenvalues 
are complex, i.e. Im(Xi) 7^ 0. 

Lemma 2.2. The solution of 2-D REP remains smooth for all time if eigenvalues are 
initially complex. Moreover there is a global invariant given by 

d 2 - Bo 2 

(2.8) ^- + 2klnp = Const. 

P 
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Proof. To obtain the desired global invariant we set q := d 2 , then from the equations 
we deduce 

dq d' q 
-^- = 2d- = -2k + (3p+-. 
dp p p 



Integration gives 



Q 
P 



— (3p + 2klnp = Const, 



which leads to ( j2.8|) . The boundedness of d follows at once, since for negative /3's, 

d 2 < max {Const. p — 2klnp + /9p 2 } =: C\. 
In particular, substitution of the lower bound, d > —C\ into the mass equation gives 



P < C lP 

which yields the desired upper-bound for the density, p(-,t) < p (a)e Clt . 



□ 



Remark. More precise information about the large time behavior is available from phase 
plane analysis. According to (p^7|), the zero level set of d' = is the hyperbola Q : = 
kp — (d 2 + /3p 2 )/2 = 0, with a right branch passing critical point (0, 0) and a left branch 
is located in the left half plane, p < 0, see Figure \2.1[ 



Q=0 




\ 




Figure 2.1. Zero level set of d' = 0. Complex spectral gap 



The trajectory on the plane p < does not affect the solution behavior in the region 
p > since p = is an invariant set, governed by 

d 2 .. . dn 
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Note that (0,0) is the only critical point of the autonomous ODE system ( |2.6| ), (2.7) 
on the right half phase plane, and that the vector field in {(p, d), Q < 0,d > 0} is 
converging to the critical point (0,0). It follows that for global smoothness it suffices to 
control the divergence d from below in the region {(p,d), Q < 0,d < 0, p > 0}, and to 
control the density from above in the region {(p, d), Q > 0}. 
For the former case we have, recalling that (3 < 0, 



p) - + 2p 2 - ' 



and its integration along particle path gives 



d>[kt+j-)p. 



This shows that the divergence d is bounded from below, and in particular, it becomes 
positive for large time since the density is positive. To the upper-bound for p in the region 
Q > 0, where d(t) > d (a), we substitute this estimate into the mass equation yielding 



P < -pd {a). 

This clearly gives the upper-bound for the density p < p (a)e~ dot . 



2.2. Real Spectral Gap. When (3 > the initial spectral gap is real, and there are two 
cases to be considered: 

Subcase 1: (3 = when the eigenvalues are equal, i.e. Ai(0) = A 2 (0). 

In this case the zero level set d' = becomes a parabola passing through the only 



critical point (0,0), see Figure p. 2], and one can repeat similar augments to our phase 
plane analysis in the previous case of distinct real roots. Note that the global invariant 
(P.8I) becomes 



d 2 

h 2klnp = Const. 

P 
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Figure 2.2. Zero level set d! = 0. Real spectral gap. 



Subcase 2: (3 > when the eigenvalues are initially real. 



Lemma 2.3. // eigenvalues of VUo are real, then the solution of 2-D REP remains 
smooth for all time if and only if 



Ai(0) + A 2 (0) >g(p ), 



where 



g(p):=sgn(p-j)Jpx (f( P ) - F fj\ ), F{p) = (3p - 2klnp. 



Proof. The system (f2~l f )- (|2.7l) has two critical points on the phase plane: 0(0,0) and 
A(f ,0), see Figure |^. 
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L(p*,d*) 



The coefficient matrix of the linearized system around (p*,d*) is 

-d* -p* 
k-pp* -d* 

A simple calculation gives the eigenvalues of L, 

A± = -d* ± y/p*(Pp*-k). 

At (0, 0), we have Ai = A2 = and hence (0, 0) is a non hyperbolic critical point. Another 

critical point, A(^-,0), is a saddle since Ai j2 = ±<J~^k. We shall use the above facts to 

construct the critical threshold via the phase plane analysis. 

Assume the seperatrix enters (leaves) A along the line d = s(p — 2^). Upon substitution 
into the linearized system around A, i.e., 



P 



-jd, d 



one can obtain 



2 



Thus, two seperatrixes leave/enter A along the directions 

6\ = — arctg and 6 2 = arctg ^ ^ 



In the phase plane the zero level set of d' = is an ellipse, see Figure |2~3| , 
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Let 7 S (A) be the portion of the stable manifold of the system coming into A from {d < 0}. 
In order to prove the existence of a critical threshold it suffices to show that 7 S (^4) can 

only come from O. Let B be the lowest point of the ellipse with coordinates y^, 

and PB be a horizontal line intersecting with p = at ^0, — J . According to the vector 

field inside the ellipse we see that the trajectory ^ S {A) can only come from the area OPB 
by crossing the curve OB. Note that the vector field on PB is going outside OPB and 
p = is invariant. Thus all trajectories in the area OPB originate from O. Therefore 
7 S (^4) can only originate from O (as t — > — oo) and becomes a portion of one unstable 
manifold of O. By symmetry we can show that the unstable manifold of the system issued 
from A entering {d > 0} will end through the portion {d > 0} at O. 
Thus the critical curve <7 : IR + — ^ IR is the one defined as 

{CM), d = g(p)} = ls (A)- 
In order to have a precise formula for g we need to use the global invariant ([2.8|), i.e., 

h 2klnp = Const. 

P 

Thus all trajectories can be expressed as 

d 2 = p[C(a) + F(p)} 

with 

C{pt) = ^ - F(po), F(p) := Pp - 2k\np. 
Po 

Note that F(p) is a convex function and min p> o F = F{^j) = 2k[l — ln{2k/(3)]. Due to 
the symmetry the homoclinic connection is possible when the trajectory pass (po, 0) with 
n - < ^ and converging to (0, 0) as t — > ±oo, i.e. the initial data must satisfy < p < 2 j 



and 

The seperatrices passing through (^,0) correspond to C(a) = —F(M-). The stable man- 
ifold 7s (A) can be written as d = g(p) for < p < oo, where 



g(j>)=Bga[p-?P\ Jp(f{p)-f( 2 4) ). 



PJ V V \PJ 

It remains to prove that the initial data satisfying do < g(po) always lead to finite time 
breakdown. 

First, in the region |(p, d), d < —Jp(F(p) — there must exist a finite time 

T x > such that p{T x ) > f for p < f ( take T x = for p > f ) since p' > 0. Therefore 
p{t) > p(Ti) for t > Ti and 



p> = - p d > Px jp[ p(T x ) - F [ j ) ) for t > Tx. 
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Integration over [Ti,t] gives 



yW)> ; V/ ^' t> Tl . 



Thus the solution must becomes unbounded before the time 

9 



yJpiTMFWTtf-Fif)) 
Second, we consider the trajectories in the region 



(P,d), P>j, \d\ < (F(p) - F 

Note that at finite time the trajectory must enter the subregion {(p,d), d < 0} through 
the left point (p*,0) identified as 

2k 

d 2 = p[F(p) - F(p*)], P>P*>j- 

This combined with the Riccati-type inequality 

df < -d 2 /2 

ensures the breakdown at finite time. This completes the confirmation of the curve d = 
g(p) as a critical threshold. □ 

Proof of Theorem 1.1: It suffices to summarize the above cases with 

Po 

being taken into account. Clearly the case (3 < and f3 — correspond to the set 

{(po, M ), r <o}. 

For (5 > 0, i.e., T > we rewrite the critical threshold as 



2k / / /2k 
do = sgn(p -JfNPo ( F (Po) ~ F [-J 



= sgn(r - 2kp )^T - 2/cpo + 2kp ln 1 
where we have used the relation F(p) = j3p — 2klnp and 

F(po) = 2klnp , 

Po 

KD— »(¥)■ 

This completes the proof of Theorem 1.1. □ 
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3. 2D REP WITH NON-ZERO BACKGROUND 

This section is devoted to the study of the REP with nonzero background c > 0, for 
which the velocity gradient tensor, M = VU, solves 

(3.1) ^M + M-^-c], 

(3.2) j^p + ptrM = 0. 



Again, using the Spectral Dynamics Lemma presented in [l7j the spectral dynamics of M 
is governed by 

k 

A2 + A2 = 2^ _c )' 
p' + p{\ 1 + \ 2 ) =0. 

As in the zero background case ,the difference rj := X 2 — Ai is proportional to the density 
along the particle path in the sense that 

77(f) 770(a) 2 

, . = 7 — r, at IK . 

p(t) po{a) 
Further manipulation gives a closed system 

(3.3) p' = -pd 

d 2 + Bo 2 

(3.4) d' = k(p-c)-—^=:Q, 0:=t>g/pjj. 



Once again the dynamics of (pT3|), ( |3.4| ) is influenced by the choice of 3. We proceed to 
discuss the solution behavior of (|3.3| ), ( p.4| ) by distinguishing two cases: 

(1) (5 < 0, the spectral gap is complex; 

(2) 8 > 0, the spectral gap is real. 



3.1. Complex Spectral Gap. We first discuss the case 3 < 0, which corresponds to 
the case that the eigenvalues are initially complex. 

Lemma 3.1. Assume that the eigenvalues are initially complex with im(Aj(0)) 7^ 0. Then 
the solution of ( \c>.Z\ ), ( \3-4\ ) remains smooth for all time. Moreover, there is a global 
invariant in time, given by 



V(p, d) = p- 1 d 2 - 3p 2 + 2kpln (0 + 2ck 



Proof. A straightforward computation yields V = along the 2D REP solutions, which 
implies that the curves V = Const are invariants of the flow. As before, for negative 3's 
we have 

d 2 < max {Const.p - 2kpln (j^j - 2ck + /5p 2 | < C 2 
and the bounds of d (and hence of p) follow. □ 
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Remark. 

1. As before, a more detailed information is available in this case by a phase plane 
analysis. If the eigenvalues are initially complex, then one has /3 = T]q/Pq < 0. The zero 
level set of d! 



Q = becomes a hyperbola, see Figure [3J. 




Figure 3.1. Zero level set d! = 0. Complex spectral gap with nonzero background. 



The intersection of its right branch with d — is the rest point A = (p*,0) of the 
system, where 

, _ k [k 2 2ck 
P -(3 + \j^-lf 
The coefficient matrix of the linearization around (p*,0) is 

w» = (*_V f)- 

Its eigenvalues satisfy 

A 2 = p*{(3p* —k) = -p*^/k 2 -2ckp < 0. 

Hence such critical point is a non hyperbolic equilibrium. The nonlinear effect plays 
essential roles in the solution behavior. In order to locate the possible critical threshold, 
we first study the solution around (p*, 0). Set n = p — p*, we then have 

(3.5) n = —p*d — nd, 

(3.6) d! = yjk 2 - 2ck(3n - y - |n 2 . 

It is easy to see that the flow governed by the linear part stays on the ellipse 

^k 2 - 2ckf3n 2 + p*d 2 = Const. 
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In order to capture the dynamics of the nonlinear system in the neighborhood of the 
critical point (n, d) = (0, 0), we employ the polar coordinates of the form 

rcosO 



n 



d 



{k 2 - 2ck(3y/ 4 '' 
—rsinO 



P 



Careful calculation with these polar coordinates, yields that (|3.5| (- (|3.6|) recast into the 
form 



(3.7) 
(3.8) 
where 



r' = R{r,B), 

9' = -^{k 2 - 2ck[3) 1 / 4 + 0(r, 9), 



R(r, 9) 
Q{r,9) 



r 2 sin9 

~ 1 WW 



1 + 



kcos 2 9 



^k 2 - 2ck(3 



rcos9 



2^^/k 2 - 2ck(3 



^k 2 - 2ckf3sin 2 9 - (3p*cos 2 9 



When r is sufficiently small 9' is strictly negative. The pleasant implication of this is that 
the orbits of system ( |3.5| ), ( |3.6| ) spiral monotonically in 9 around the (p*, 0). But the even 
power of r 2 does not tells the stability property of the critical point. 

Observe that if (n(t),d(t)) is a solution, so is (n(—t), —d(—t)). Such symmetry implies 
that there is a center in the neighborhood of (p*, 0). 

In order to clarify the global behavior of the flow around such center, we appeal to the 
global invariant 

V(p, d) = p" 1 d 2 - (3p 2 + 2kpln + 2ck 

we claim that V is positive definite which serves as a (majoration of) Lyapunov functional, 
to this end, we consider the function H(p) := —f3p 2 + 2kpln (^) + 2ck, which is convex 
and takes its minimum at p m i n , satisfying 



In 



Pu 



2c 



1 + k P » 



Observe that since /3 < the function 

= 0, which when combined with the fact 



is an increasing function in p > and h(p 
that h(p*) > verifies that 



< p r 
Indeed, for p* we have 



P 



P~ 1 [k-y/k 2 -2ck/3]. 



Hp* 



1 



-p + In 

k 



= y/l - 2ck~ l p -ln(l+ v 7 ! - 2cfc- 1 /?) > 0. 
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H(Pmin) = PP^in ~ 2kp min + 2ck = (3(p mhi - p*){p m \n ~ P*) > 0, 

where p* = p~ l [k + ^/k 2 - 2ckft\. 

The invariant curves, V = Const., represent of course, the bounded periodic orbits 
containing (p*, 0). 

3.2. Real Spectral Gap. We divide the region ft G [0, oo) into subregions depending 
on the number of critical points on the phase plane, and then study the solution behavior 
with ft in each sub-region. The solution behavior depends strongly on the number of 
critical points and their stability property. 

Let (p*, d*) be a critical point of the system, then the coefficient matrix of the lineariza- 
tion around (p*,d*) reads 

W) = { k -%* :£)■ 

Its eigenvalues are given by 

(3.9) A = -d* ± sjp*{ftp*-k). 

We now discuss subcases distinguished by the number and type of the critical points 
as P changes. 

• ft — 0, then the zero level set d' = Q = is a parabola, d 2 = 2k(p — c), intersecting 
with d — at (p*, d*) — (c, 0). From ( |3.9| ) we see that at this point the eigenvalues of 
L are A = ±\/cki, pure imaginary number, the critical point (c, 0) is non-hyperbolic. 
The stability property of this critical point has to be determined by taking into 
account of the nonlinear effect. 

• < ft < y c - The zero level set Q = is an ellipse, located on the right half plane 
p > 0. There are two critical points (p*,d*) = (p*,0) with 

„ k [k 2 2k~c 

The associated eigenvalues of L are 

A(p^) = ±yW^ 2 -2c^, A(p* 2 ) = ±^/ P Wk 2 -2ck^. 

Therefore (pi, 0) is a center of the linearized system and (pSj, 0) is a saddle, see Figure 
3~2|. 




p 



Figure 3.2. Critical points in p — d plane. Real spectral gap with nonzero background. 



Possible bifurcation as (3 changes from to ^- may be responsible for the compli- 
cated solution structure in this regime. 
/? = 7T-. The zero level set Q = 0, i.e., 



d 2 + P( P -^) 2 



degenerates to a single point (p*,d* 
eigenvalues. 
(3 > In this case 



'§,0), the only critical point with zero 



2c ' 



Q 



d 2 + (3 ( p - - ) + 2kc - 



k 2 , 
< — - kc < 0. 

" 2(3 



There is no critical point at all in the finite phase plane. 
The solution behavior distinguished by above cases is given in the following lemmata. 

Lemma 3.2. If Ai(0) = A2(0). Then the solution of ( \3.3j) , (3.4 ) remains smooth for all 
time, indicated by the global invariant 

d 2 + 2ck 



(3.10) 



+ 2fclnp = Const. 



Proof. The assumption amounts to (3 = 0. As discussed above, (c, 0) is the only critical 
point and the center of the linearized system. In order to find the global invariant we set 
:) 2 + cd 2 . Along the trajectory ^ = U(x,t) 



R(t) := k(p - 
(3.11) 

/From the mass equation it follows that 



— R(t) = 2k(p - c)p + 2cdd' = -d[k(p - c) 2 + R(t)]. 



d = -'- 
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which when inserted into the relation ( p.ll| ) gives 

<LR k(p - cf R 

~T = - + -■ 

dp p p 

Integration gives 

— + ^— + 2cklnp — kp = Const, 
P P 



which leads us to the global invariant as asserted in (|3.10| ). This global invariant is 
compact and ensures that both divergence d and the density p remains bounded as time 
evolves. □ 

We leave the case < (3 < later and study the critical case (3 = j-- 

Lemma 3.3. If A 2 (0) — Ai(0) = y^Po- Then the solution of ( \3.3j) , ([3.4 ) always develops 
finite time breakdown unless the initial data lies in the set 

d 2 4- 2ck k 

(p, d) G IR + x IR + , p + 2klnp = 2kln(2c) 

p 2c 

Proof. The given assumption is equivalent to the case P = In this case the divergence 
always decrease except for at the critical point (2c, 0) since 

* = -£-*("- *)'<„, 

2 4c ~ 

In order to clarify the solution behavior we proceed to obtain the global invariant. Set 
q := d 2 , one then has 

dq _ 2dd' _q + Pp 2 - 2k(p - c) 
dp p' p 
Solving this equation we obtain 

- = Pp — 2klnp h Const. 

P P 

Therefore we come up with a global invariant 

d 2 + 2rk 

(3.12) — pp + 2klnp = Const. 

P 

The only trajectory converging to the critical point is realized by a half trajectory con- 
verging to (2c, 0) from the first quadrant. For all other trajectories not passing the critical 
point (2c, 0), the rate d' is strictly negative. The divergence will become negative at finite 
time even if it is initially positive, which when combined with the Riccati-type inequality 
d' < —d 2 /2 confirms the finite time breakdown. □ 



We now look at the case P > 



2c' 



Lemma 3.4. Assume that the eigenvalues are initially real and | A2 (0) — Ai (0) | > y ^po(a). 
Then the solution of ( \3. $) , ($-4) always develops finite time breakdown. 
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Proof. The given assumption is nothing but the inequality /3 > Note that there is 
no critical point in the finite phase plane, actually Q remains negative for all time. The 
solution must develops breakdown in finite time. In fact from 

d 2 (3 ( k\ 2 kc ( k 



we find that 

d'<-8 with 6:=j(p-£)>0. 

This ensures that d must become negative beyond a finite time T , say T > maxjy, 0}. 
The d— equation ( |3.13|) also gives 

d! < -d 2 /2, 

whose integration over [T ,t] leads to 



l-\d(T a )(t-T ) 

nro /.. — 

d(T )- 



Hence the solution must break down at a finite time before T — -n^. □ 



Finally we conclude this subsection by discussing the delicate case < (3 < j-. Set 

G(p, p\(3) := (3( P - p*) - 2kln (£\_ 2 _±+ 2 _± 

\P J P P 



with p* = (3- l [k ± ^k 2 - 2ck(3] being the p-coordinate of the intersection point of the 
trajectory with p axis. 



Lemma 3.5. Assume that the real eigenvalues satisfy < |A 2 (0) — Ai(0)| < y-^p (a). 

Then for any (3 E (0, the solutions of \3. 3j ), ([3.4) remains smooth for all time if and 
only if 

|Ai(0) + A 2 (0)| < ^p G(po,p* 2 ,(3o) for p < p\ 



and 



Ai(0) + A 2 (0) = \/poG(po,p2*,Po) for p > p*. 



Proof. The assumption tells us that /3 < j-. In this case there are two critical points in 
the phase plane A = (p*,0) and B = (pj^O), see Figure |3T2| . B is a saddle, whose two 
manifolds pass enclosing the critical point A, which is a center of the linearized system. 
Let W g (B) denote the stable manifold coming from the region {(p,d), p < p* 2) d < 0} and 
W U (B) the unstable manifold entering into {(p,d),p < p* 2) d> 0}. To prove the results 
stated in the theorem it suffices to show that for any (3 G (0, ^) such that 

w u (B)nw s (B) 

is not empty, i.e. there always exists a saddle connection (homoclinic orbit). 
This follows from the continuity argument supported by the following facts: 
(1) Both W U (B) and W S (B) pass through the segment OA with flow going downward 

since d! < and p' = on OA, see Figure |3~2 

The level curve d' = is an ellipse with upper vortex P located at (4, 0). Let P\ denote 
the intersection of tangent line of the ellipse Q = through P with the axis p = 0. The 
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vector field inside the ellipse shows that W U (B) must escape the ellipse from the curve 
PA. Note that the trajectories on PP 1 and AP enters into the region PAOP\, and the 
axis p = is an invariant set. These facts ensure that W U (B) must enter the region d < 
through OA. Similarly we can show that W S (B) for p < p\ must enter the region d > 
through OA. 

(2) As (3 increase in (0, the point W U (B) fl OA moves to the right and the point 
W S (B) n OA moves to the left. 

We prove the claim for W U (B) fl OA, and the case for W S (B) n OA follows similarly. 
The claim follows from the following two observations: 

(i) The slope of the unstable manifold W u {B,f$) at (pl,0) is d p d\ p=p * = \-(p* 2 ,f3) and 
the eigenvalue \_(pl,(3) is increasing in (3; Indeed 

dp k\-{p* 2 ,(3) I y ^k 2 - 2ck(3 J 

(ii) W u (B,l3i) does not intersect with W u (B,f3 2 ) for /3i ^ /3 2 ; 

As done previously we can find the global invariant of the system 



2ck 

f3 p — 2klnp h Const 

P 



d 2 = p 

from which the left branch of the unstable manifold of B can be explicitely expressed as 

< P < P* 2 . 



A careful calculation gives 



(3(p - p*) - 2kln(p/p*) - — + ^ 

P P2 



dd p . 



which ensures the claim (ii). 

(3) Let p u {(3) be the p coordinate of the point W U (B, (3)C\OA and p s (/3) be p-coordinate 
of the point W S (B, (3) n OA. We claim 

lim p u ((3) < lim p s {(3). 

In fact from the expression of seperatrices 

d 2 = pG(p,p* 2 ,(3), p<p*, 
we find that the p-coordinate of points W U / S (B) fl OA satisfy 

G(p,pl,P) = 0. 

Note that 

dG (3 dG 

9j = f(p-Pi)(p-P*), ^=P~P2- 



Thus we have for < p < p\ 



dp dG .dG p 2 ,, . n 
Oi3 = 'di3 / ^ = J iPl - p)> °- 



This confirms the above assertion. 
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Combining the above observations we conclude that there exists a (3 G (0, ^) for 
which a saddle connection exists. It remains to show that as (3 changes in the region 
(0, y) the above saddle connection is preserved. Observe that if (p(t),d(t)) is a solution, 
so is (p(— t), —d(—t)). Such symmetry prevents the occurrence of the possible bifurcation 
when (3 changes. 

Using the nonlinear terms in the equation and the vector field we can show for the 
initial data outside the closed curve — saddle connection — the solution always develops 
finite time breakdown, details are omitted. □ 



Proof of Theorem 1.2: Summarizing the results stated in the above lemmas we see 
that the case (3 < and (3 = corresponds to the set 



Si = i (po, d , T ) 



r n < and 



d > if po = 

d arbitrary if p > 



since T = (3p^. The case < (3 < ^ corresponds to < T < ^Po^ ^ ne divergence is 
required to satisfy the critical threshold condition 

|do| < VpoG(po,P* 2 ,P), < po < p* 2 
and d = a/ PoG(p , p* 2 , (3) for p > p\. Using T = (3p\ and 

2ck 2cp 



p* 2 = (3~ l [k + ^k 2 -2ck(3] 



V - y/W^2ckf3 PQ - s JpJ^: 



one has 



p G{po, p* 2 , (3) = po 



(3(p - pi) - 2kln 



= r i 



2c 



Po - \/Po~ — 




2ck 2ck 
Po P* 2 



2kp ln 



In 2 2cT 

Po ~ \/P o - — 
2c 



2ck + k \ po 



2 2 C r 

Po 



k 



= To-2ck-2k\ pl- 



2cr 
k 



— 2kpoln 



In 2 2cr » 

Po - \/Po — 
2c 



which leads to the critical threshold described by the set S2. The set S3 can be figured 
out in a similar manner. 

□ 

ACKNOWLEDGMENTS 



Research was supported in part by ONR Grant No. N00014-91-J-1076 (ET) and by 
NSF grant #DMS01-07917 (ET, HL). 



Critical Thresholds in restricted Euler-Poisson equations 



23 



References 

[1] Olu§ Boratav and Richard Pelz, On the local topology evolution of a high- symmetry flow, Phys. 
Fluids 7 (1995), 895-897 and 1712-1731. 

[2] U. Brauer, A. Rendal, AND O. Reula, The cosmic no-hair theorem and the non-linear stability 
of homogeneous Newtonian cosmological models, Class. Quantum Grav. 11 (1994), 2283-2296. 

[3] Brain J. Cantwell, Exact solution of a restricted Euler equation for the velocity gradient tensor, 
Phys. Fluids A, 4 (1992), 782-793. 

[4] M. Chertkov, A. Pumir AND B. Shraiman, Lagrangian tetrad dynamics and phenomenology of 
turbulence, Phys. Fluids A, 11 (1999), 2394-2410. 

[5] G.-Q. Chen and D. Wang, Convergence of shock capturing scheme for the compressible Euler- 
Poisson equations, Comm. Math. Phys. 179 (1996), 333-364. 

[6] J. Dolbeault and G. Rein, Time- dependent rescaling and Lyapunov functionals for the Vlasov- 
Poisson and Euler-Poisson systems, and for related models of kinetic equations, fluid dynamics and 
quantum physics, Mathematical Models and methods in Applied Sciences. 

[7] S. Engelberg, Formation of singularities in the Euler-Poisson equations, Physica D 98 (1996), 
67-74. 

[8] S. Engelberg, H. Liu and E. Tadmor, Critical Thresholds in Euler-Poisson Equations , Indiana 

Univ. Math. J., 50, (2001), 109-157. 
[9] P. Gamblin, Solution rgulire temps petit pour I' quation d 'Euler-Poisson, Comm. Partial Differential 

Equations 18 (1993), 731-745. 
[10] I. Gasser, C.-K Lin and P. A. Markowich, A review of dispersive limits of (non)linear 

Schrddinger-type equations, preprint, 2000. 
[11] Y. Guo, Smooth irrotational flows in the large to the Euler-Poisson system in IR 3+1 , Comm. math. 

Phys. 195 (1998), 249-265. 
[12] D. Holm, S. F. Johnson, and K.E. Lonngren, Expansion of a cold ion cloud, Appl. Phys. Lett. 38 

(1981), 519-521. 

[13] J.D. Jackson, Classical Electrodynamics, 2nd ed., Wiley, New York, 1975. 

[14] S. Junca and M. Rascle, Relaxation of the isothermal Euler-Poisson system to the drift- diffusion 

equations, Quart. Appl. Math. 58 (2000), 511-521. 
[15] T. Luo, R. Natalini, and Z. Xin, Large time behavior of the solutions to a hydrodynamic model 

for semiconductors, SIAM J. Appl. Math. 59 (1999), 810-830. 
[16] H. Liu and E. Tadmor, Critical thresholds in a convolution model for nonlinear conservation laws, 

SIAM J. Math. Anal. 33 (2001), 930-945. 
[17] H. Liu and E. Tadmor, Spectral dynamics of velocity gradient field in restricted flows, to appear 

in Comm. Math. Phys. 

[18] T. Makino, On a local existence theorem for the evolution of gaseous stars, Patterns and Waves 
(eds. T. Nishida, M. Mimura and H. Fujii), North-Holland/Kinokuniya, 1986, 459-479. 

[19] P. A. Markowich, A non-isentropic Euler-Poisson Model for a Collisionless Plasma, Math. Meth- 
ods Appl. Sci. 16 (1993), 409-442. 

[20] T. Makino, and B. Perthame, Sur les solutions a symetrie spherique de I'equation d'Euler- 
Poisson pour I'etoiles gazeuses, Japan J. Appl. math., 7 (1990), 165-170. 

[21] P. Marcati and R. Natalini, Weak solutions to a hydrodynamic model for semiconductors and 
relaxation to the drift- diffusion equation, Arch. Rat. Mech. Anal. 129 (1995), 129-145. 

[22] P. A. Markowich, C. Ringhofer, and C. Schmeiser, Semiconductor Equations, Springer, 
Berlin, Heidelberg, New York, 1990. 

[23] T. Makino and S. Ukai , Sur I'existence des solutions locales de Equation d 'Euler-Poisson pour 
I'volution d'toiles gazeuses, J. Math. Kyoto Univ. 27 (1987), 387-399. 

[24] B. Perthame, Nonexistence of global solutions to the Euler-Poisson equations for repulsive forces, 
Japan J. Appl. Math. 7 (1990), 363-367. 

[25] J.B. Thoo and J.K. Hunter, Nonlinear wave propagation in a one- dimensional random medium, 
preprint (2001). 



24 



H. Liu and E. Tadmor 



[26] P. Vieillefosse, Local interaction between vorticity and shear in a perfect incompressible flow, J. 

Phys. (Paris) 43, (1982), 837. 
[27] D. Wang, Global solutions and relaxation limits of Euler-Poisson equations, to appear in Z. Angcw. 

Math. Phys. 

[28] D. Wang, Global solutions to the equations of viscous gas flows, to appear in Pro. Royal Soc. 
Edinburgh: Section A. 

[29] D. Wang Global solutions to the Euler-Poisson equations of two- carrier types in one- dimension, Z. 

Angew. Math. Phys. 48 (1997), 680-693. 
[30] D. Wang, G.-Q. Chen, Formation of singularities in compressible Euler-Poisson fluids with heat 

diffusion and damping relaxation, J. Diff. Eqs. 144 (1998), 44-65. 

UCLA, Mathematics Department, Los Angeles, CA 90095-1555. 
E-mail address: hliu@math.ucla.edu 

UCLA, Mathematics Department, Los Angeles, CA 90095-1555. 
E-mail address: tadmor@math.ucla.edu 



